RANDOM REGULARIZATION OF BROWN SPECTRAL 

MEASURE 

b 

^ ; PIOTR SXTADY 

> 

O' Abstract. We generalize a recent result of Haagerup; namely we 

, show that a convolution with a standard Gaussian random ma- 

trix regularizes the behavior of Fuglede-Kadison determinant and 
, Brown spectral distribution measure. In this way it is possible to 

establish a connection between the limit eigenvalues distributions 
of a wide class of random matrices and the Brown measure of the 
corresponding limits. 
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1. Introduction 

>■ \ The problem of determining the joint distribution of the eigenvalues 

of a random matrix G A^Ar(£°° _ (fi)) with a given distribution 
of entries is usually very hard and has explicit solutions only for very 
limited cases. A partial solution of this problem is to consider the limit 
iV — > oo and to hope that it is an easier problem than the original one. 
O ' Questions of this type can often be solved by Voiculescu's theory 

c-| ■ of free probability [VDN|: if a sequence (A( N ') of random matrices 



X 



converges in ^-moments to some x, where x is an element of a non- 
commutative probability space (A, 0), then some properties of matrices 
A( N ) ( e .g. independence of some entries) imply useful properties of the 
limit object x (e.g. some kind of freeness). If random matrices A^' are 
normal then the empirical distribution of their eigenvalues converges to 
the spectral measure of x. This approach turned out to be successful in 
determining the limit eigenvalues distribution in many cases (see e.g. 

The situation is much more complicated if random matrices and 
the limit object x are not normal. In this case the spectral measure 
of x has to be replaced by a more complicated object, namely by the 
Brown measure of x. The problem of determining the Brown measure 
of a given operator is still difficult, but it is much easier than the orig- 



inal question about the eigenvalues of a random matrix [[Larj , |HL| , |BL 
However, since the Brown measure does not behave in a continuous 
way with respect to the topology given by ^-moments, the distribution 
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of eigenvalues of A^ N ' does not always converge to the Brown measure 
of x. 

Surprisingly, in many known cases when we consider a "reasonable" 
or "generic" sequence of random matrices the sequence of their Brown 
measures converges to the Brown measure of the limit (cf ||BL|| ) . In 
these cases, however, the convergence was proved in this way, that the 
distribution of eigenvalues of was calculated by ad hoc methods 
and nearly by an accident it turned out to converge to the Brown 
measure of x. Therefore one of the most interesting problems in the 
theory of random matrices is to relate the asymptotic distribution of 
eigenvalues of a sequence of random matrices with the Brown measure 
of the limit. 

In this article we show that it is possible to add a small random 
correction to a sequence of random matrices (A^ N ') which converges in 
^-moments almost surely to some element x G (A, 0) in such a way that 
the new corrected sequence still converges to x and that furthermore 
the empirical eigenvalues distributions of the new sequence converge 
to the Brown measure of x almost surely. A slightly different version 
of this result was proved recently by Haagerup |[Haa|| and plays a key 
role in his proof of existence of invariant subspaces for a large class 
of operators. The random correction used by him is a matrix Cauchy 
distribution, the first moment of which is unbounded, what makes it 
unpleasant for applications. For this reason Haagerup's estimates of 
the correction were in the C p norm with < p < 1. The random 
correction considered in this article has a nicer form of a Gaussian 
random matrix and for this reason we are able to find better estimates 
for the correction, namely in the operator norm. 

The Gaussianity of the correction allows us also to find the limit 
empirical eigenvalues distribution of a wide class of random matrices, 
which include both the well-known examples of the matrix G^ N ' with 
suitably normalized independent Gaussian entries (the limit eigenval- 
ues distribution was computed by Ginibre ||Gin|| in the sense of density 
of eigenvalues, the almost sure convergence of empirical distributions 
was proved in unpublished notes of Silverstein and later in more gen- 
erality by Bai ||Bai|| ; the Brown measure of the limit was computed by 
Haagerup and Larsen | |Lar| , [Hh|| ) , the so-called elliptic ensemble (the 
limit eigenvalues distribution was computed by Petz and Hiai ||PH|| and 
the Brown measure of the limit was computed by Haagerup and Larsen 
[|Lar| , |HL|1 ) and new examples for which the eigenvalues distribution was 
not known before and which are of the form G^ N ' + A^ N \ where entries 
of G^ and A^ are independent (the Brown measure of the limit of 
such matrices was computed by Biane and Lehner ||BL|1 ). 



RANDOM REGULARIZATION OF BROWN SPECTRAL MEASURE 



3 



Results of this article can be also applied ||Siii|| to show that DT 
operators (which were introduced recently by Dykema and Haagerup 
[PH|| ) maximize microstate free entropy [|Vo2|| among all operators hav- 



ing fixed Brown measure and the second moment. 

Our method bases on the observation that if a sequence of random 
matrices converges in ^-moments to x then the Fuglede-Kadison 
determinants A(A^) converge to A(x) as well if we are able to find 
some bottom bounds for the smallest singular values of A^ N \ Since 
the random correction considered in this article is given by a certain 
matrix-valued Brownian motion, hence we are able to write a system of 
stochastic differential equations fulfilled by the singular values. Unfor- 
tunately, finding an exact analytic solution to a non-linear stochastic 
differential equation is very difficult. We deal with this problem by 
proving a certain mononicity property of our equations and hence we 
are able to find appropriate bottom estimates for the singular values. 

2. Preliminaries 

2.1. Non— commutative probability spaces. A non commutative 
probability space is a pair (.4, 0) , where A is a C*-algebra and (p is 
a normal, faithful, tracial state on A. Elements of A will be referred 
to as non-commutative random variables and state (J) as expectation 
value. The distribution of x G A is the collection of all its ^-moments 
{(f)(x si ■ ■ ■ x Sn )), where Si, . . . , s n G {1, *}. 

2.2. Fuglede-Kadison determinant. Let a non-commutative prob- 
ability space (.4, 4>) be given. For x & Awe define its Fuglede-Kadison 
determinant A(x) by (cf |[FK|| ) 

A(x) = exp [0(ln |x|)] . 

2.3. Brown measure. Let a non-commutative probability space 
(A, 0) be given. For x G A we define its Brown measure ||Bro | to 
be the Schwartz distribution on C given by 

1 / d 2 d 2 \ 

One can show that in fact fi x is a positive probability measure on C. 

Example. The Brown measure of a normal operator has a particularly 
easy form; let x G A be a normal operator and let E denote its spectral 
measure: 



x 



/ z dE(z). 
Jc 
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Then the Brown measure of x is given by 

fi x (X) = <f>[E(X)\ 
for every Borel set X C C and the following holds: 

(t>[x k {x*) 1 ] = [ z k z l d(i x (z). 
Jc 

2.4. Random matrices. We have that (.Mjv,trjv) is a non- 
commutative probability space, where A4n denotes the set of all 
complex- valued N x N matrices and trjy (which for simplicity will 
be also denoted by tr) is the normalized trace on A4 N given by 

tr N A = ^-TtA iorAeM N , 
N 

and Tr denotes the standard trace. 

The below simple example shows that for finite matrices the Fuglede- 
Kadison determinant A and the usual determinant det are closely re- 
lated and gives heuristical arguments that for every Borel set X C C 
the Brown measure /jl x (X) provides information on the joint "dimen- 
sion" of "eigenspaces" corresponding to A G X. 

Proposition 1. The Fuglede-Kadison determinant of a matrix A G 
.Mat with respect to a normalized trace tr is given by 

A(A) = y\detA\. 

The Brown measure of a matrix A G Mn with respect to the state 
tr is a probability counting measure 

1 N 

i=i 

where X 1 , . . . , \ N are the eigenvalues of A counted with multiples. 

In the following we will be interested in studying the random measure 
uj i— > HAioj) for a random matrix A G A4n(£°°~ (^)) • This random 
measure is called the empirical distribution of eigenvalues. 

We will use the following convention: we say that random matrices 
A, B G M.n{£°°~{Q)) are independent if the family of entries of A and 
the family of entries of B are independent. 

2.5. Convergence of moments. Let a sequence 
•MAr(£°°~(f2)) of random matrices, a non-commutative probability 
space (A, 4>) and x G A be given. We say that the sequence A^ 
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converges to x in ★-moments almost surely if for every n G N and 
Si,...,s n G {1,*} we have that 

lim ti N [(AW)" 1 • • • = 0(a; Sl • • ■ x s ") 

holds almost surely. 

Let a sequence G Mn(£°°~(VI)) of random matrices, a non- 

commutative probability space (A, 4>) and x G ^4 be given. We say that 
the sequence converges to x in expected ★-moments if for every 

n G N and s\, . . . , s n G {1, ★} we have that 

lim EtT N \(AW) ai ...(AWy»] =<f>( x 'i...x 9 "). 

N~ +oo 

2.6. Discontinuity of Fuglede— Kadison determinant and 
Brown measure. One of the greatest difficulties connected with the 
Fuglede-Kadison determinant and Brown spectral distribution measure 
is that — as we shall see in the following example — these two objects do 
not behave in a continuous way with respect to the topology given by 
convergence of ★-moments. 

We say that u G A is a Haar unitary if u is unitary and 4>{u k ) = 
0((u*) k ) = for every k = 1,2,... It is not difficult to see that the 
sequence (S^) converges in ★-moments to the Haar unitary, where 
is an N x N nilpotent matrix 

•••0 0" 

1 •••0 

1 ••. : : . 

1 '•• '•• 
••• 1 _ 

Every matrix has the determinant equal to 0, while the Haar 

unitary has the Fuglede-Kadison determinant equal to 1; every matrix 
E.W has the Brown measure equal to So, while the Brown measure of 
the Haar unitary is the uniform measure on the unit circle {z G C : 

1*1 = 1}- 

The reason for the discontinuity of Fuglede-Kadison determinant 
is that the logarithm is not bounded from below on any interval [0, t}. 
However, since it is bounded from above, Fuglede-Kadison determinant 
is upper-semicontinuous. 

Lemma 2. Let AW be a sequence of random matrices which converges 
in -k-moments to a non-commutative random variable x almost surely. 
Then for every A G C 

lim sup tr In \A {N) - A| < In A(x - A) 

N^oo 



(1) 
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holds almost surely. 

Let be a sequence of random matrices which converges in ex- 
pected * -moments to a non-commutative random variable x. Then for 
every A e C we have 

limsupEtrln|A (A ° - A| < lnA(x - A). 

AT— >oo 

Proof. For each e > there exists an even polynomial Q such that 
lnr < Q(r) for every r > 

and 

ln(r 2 + e) 

Q{f) < — 2 — ^ or ever y < r < ||x||. 

Hence 

ti\n\A {N) - A| < trQ(|A (7V) - A|). 
The right-hand side converges almost surely (resp. in the expectation 
value) to (j)[Q(\x- A|)] < <f> ( ln(T ^ +e) ). By taking the limit e -» both 
parts of the lemma follow. □ 

2.7. Gaussian random matrices. We say that a random matrix 

G ( N ) = (cgO)^ e (n)) 

is a standard Gaussian random matrix if 

are independent Gaussian variables with mean zero and variance 277. 
We say that 

MW : R + - M N (£">-{n)), M^ N \t) = (M^it))^. .^ 
is a standard matrix Brownian motion if 

(»< ) ) 1 < iJ -< K .(3< > ) 1 <„< 1 v 

are independent Brownian motions which are normalized in such a way 
that the variance is given by 

t 

2N' 



2.8. Circular element. There are many concrete characterizations of 
the Voiculescu's circular element c [VPN] but we will use the following 
implicit definition. One can show that the sequence converges 
both in expected ^-moments and in ^-moments almost surely to a 
certain non-commutative random variable c II Vol! ITho . 
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2.9. Preeness. Let (A, <fi) be a non-commutative probability space 
and let (Ai)i^i be a family of unital *-subalgebras of A. We say that 
the algebras (Ai)i^i are free if 

4>(xix 2 ■ ■ -x n ) = 

holds for every n > 1, every i l5 z 2 , . . . , i n G / such that %i ^ i 2 , ii ^ 13, 
. . . , i n _i 7^ i n , and every Xi G Ai 1 , . . . ,x n G Ai n such that </>(xi) = 

■ ■■ = <p(x n ) = o(dwm)- 

Let (Xj)j e / be a family of subsets of A. We say that sets Xj are free 
if unital ★-algebras (A\g{Xi, X*}) i£ i are free. 

3. Singular values of matrix Brownian motions 

Let A G -Mjv(£°°~ (^)) be a given random matrix and be a 

standard matrix Brownian motion such that A and are indepen- 

dent. For t > we define a random matrix A t by 

^(w) = A(lo) + M {N \t,u). 

It should be understood that matrix Brownian motions (M^)jv=i,2,... 
are independent. 

If we are interested in A t for only one value of t > we can express 
A t as follows: 

A t (u)=A(io) + ViG {N \u), 
where is a standard Gaussian random matrix such that A and 

are independent. 

If x G A is a non-commutative random variable, we we can always 
extend the algebra A and find c £ A such that {x, x*} and {c, c*} are 
free and c is a circular element | VDN|| . We will denote 



X t = X + Vt c. 

Proposition 3. If sequence of random matrices \A^\ 2 converges in 
-k-moments to \x\ 2 almost surely then for every t > the sequence 
\A[ N ^\ 2 converges in -k-moments to \xt\ 2 almost surely. 

If sequence of non-random matrices \A^\ 2 converges in -k-moments 
to \x\ 2 then for every t > the sequence \A[ N ^\ 2 converges in expected 
•k-moments to \xA 
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Proof. The first part of the propositions follows under additional as- 
sumption that sup N IIA^H < 00 almost surely from recent results of 
Hiai and Petz [|HP|| . For the general case observe that since for all 
unitary matrices U,V G M. n and n G N the distributions of ran- 
dom variables ti\A^ + ^l 2 " and ti\A^ + y/tUG^V \ 2n = 
tr\V*AWu* + ^GW| 2n coincide, hence it is enough to prove the 
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first part under assumption that every matrix is almost surely 
diagonal. The method of Thorbj0rnsen can be generalized to this case 
fflTSo l. 



The second part of the proposition was proved by Voiculescu |[Vol |. 

□ 

For any t > and u G ft let Ai(t, to) > ■ ■ ■ > Ajv(i, oj) denote singular 
values of the matrix A t (uj). 

In Section |5.1| we derive stochastic differential equations for 
Ai, . . . , Xn using similar methods to those of Chan [ Cha ] and obtain 



(2) d\i(t) 



^ + £ f 1 " 237 + ? 



A? + Af 



£ m - a?) 

where £> is a standard matrix Brownian motion. 

Theorem 4. Lei and be non-random matrices of the same 
size, A^\ A^ G Mn- For n = 1, 2 Ze£ > • • • > sy be the singular 
values of the matrix A^ . Suppose that for each 1 < k < N we have 

Then for every t > there exists a probability space (O, i3, P) and 
random matrices G' 1 ', G^ 2 - 1 G .Mat snc/i that each matrix G' 1 ' 

zs a standard Gaussian random matrix (but matrices and G^ 
might be dependent) and such that 



tr f(\A^ + Vt GM(u)\) < tr f(\A^ + Vt 

holds for every u G Q and every nondecreasing function f : R — ► R. 

Proof. Let us consider a probability space (Q, B,P), a standard matrix 
Brownian motion £? : R + — > _Mjv(£°°~(fi)) and for each n G {1, 2} we 
find the solution of the system of stochastic differential equations 



(3) d\r> = &(dBii) + 



^ ■ dt d 1 iV (Aj n) ) 2 + (Af) a 

2Af> I 2iV ^ ((A W )2 _ (A (n) )2) 



together with initial conditions 

Ai n) (o, w ) = a w 

We have that for each t > and n G {1,2} the joint distribution of ran- 
dom variables A^(t), z = 1, . . . , N coincides with the joint distribution 
of singular values of the random matrix A[ n \ 
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The theorem will follow from the following stronger statement: for 
almost every uj and every t > we have 

(4) Af ] (t, u) < Xf\t, u) for every 1 < i < N. 

From Eq. (|3|) it follows that for almost every uj G Q we have that 
_ \( 2 ) h as a continuous derivative (see Section |5.1.2|) . For a fixed 



u E Q let to be the smallest t > such that @ does not hold. Trivially 
we have to > 0. There exists an index j such that A,- (to) = ^f\t 



Xj and for every % we have A^(t ) < A^(to). Eq. (S) gives us 



(5 ) ^(A» W -Af W ) =J2 ^' ? - ( ^ 



.... iv(A| - (A? - (A< 2 >)*) 



It is easy to see that if there exists at least one index 1 < k < N 
such that X^\t ) + A^ 2) (t ) then 



<0, 

t=to 



so it follows that for small d > we have Xj (t) — Aj (t) > for 
to — d < t < t . This contradicts the minimality of to- 

We define 6{t) = A, (1) (t) - Xf \t). If we replace in (§ Af } by A, (1) - S t 
then it becomes a system of non-stochastic ordinary differential equa- 
tions for 5i. If for all indexes 1 < i < N we have A^(to) = Xf\t 



o 



then 5i(to) = and the solution exists and is unique in some (backward) 
interval. This contradicts the minimality of to- □ 

Proposition 5. If E M n (C°°-(D,)) is a random matrix and X E 

C then the function W + 3 t i— > Etrln \A[ N ^ — X\ is nondecreasing. 

For x & A and X E C we have that the function M + 3 1 1— ► In A(x t ) 
zs nondecreasing and 

(6) lim In A(x t ) = In A (a?). 

Proof. We can regard .M 2iV 2 -dimensional real Euclidean space 

equipped with a scalar product (m,n) = KTrmn*. As usually we 
define the Laplacian to be V 2 = J2i<k<2N 2 ^v h i where V\, . . . , v 2 n 2 is 
the orthonormal basis of this space and D v is a derivative operator in 
direction v. 

Notice that In | det A\ = 3? In det A. We can regard det A as a holo- 
morphic function of N 2 complex variables (=entries of the matrix). On 
the other hand it is a known-fact that if f(zi, . . . , z^) is a holomorphic 
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function then the Laplacian of its logarithm is a positive measure. This 
and Ito formula imply the first part of the proposition. 

For the second part we construct a sequence (A^), where G 
A4n, such that A^ converges in ^-moments to \x\ and such that 
limjv^oo tr ln(A^) = lnA(x) and apply Lemma |2| for the sequence 
At . This shows that 1R + 3 1 1— > In A(x t ) is nondecreasing. 

On the other hand the inequality 

limsup In A(xt) < lnA(x) 
can be proved similarly as in Lemma [| □ 

4. The main result 

Theorem 6. Let A^ N > G .Mjv(£°°~(fi)) be a sequence of random ma- 
trices such that A^ converges in -k-moments to a non-commutative 
random variable x almost surely. 

For every t > we have that the sequence of empirical distributions 
i^ A W,s converges in the weak topology to fi Xt almost surely. 

There exists a sequence (t^) of positive numbers such that 

limjv^ootjv = an d the sequence of empirical distributions , 

t N y^) 

converges in the weak topology to fx x almost surely. 

Theorem 7. Let (A^) be a sequence of non-random matrices 
(A^ G AAn) which converges in -k-moments to a non-commutative 
random variable x G A, where {A, <p) is a non-commutative probability 
space. 

There exists a sequence (A^) of non-random matrices such that the 
distributions of eigenvalues /i^(jv) converge weakly to fi x and 

lim pW-iW||=0, 

where \\ ■ \\ denotes the operator norm of a matrix. 

As an illustration to the above theorems we present on Fig. p] — |] re- 
sults of a computer experiment; we plotted eigenvalues of the nilpotent 
matrix Sjy from Eq. (|l|) with a random Gaussian correction. The size 
of the matrices was iV = 100; with dashed line we marked the spectrum 
of the Haar unitary, which is the circle of radius 1 centered in 0. We 
recall that the sequence (H^) converges in ^-moments to the Haar 
unitary. As one can see if the random correction is too small then the 
eigenvalues of the corrected matrix behave like the the eigenvalues of 
and if the random correction is too big then the eigenvalues of the 
corrected matrix are dispersing on the plane. 



RANDOM REGULARIZATION OF BROWN SPECTRAL MEASURE 



11 















\ 


\ 


/ 

/ / 
/ / 




\ 


\ 

\ 

j 


\ \ 






/ 


\ \ 






/ 

/ 


\ 

\ 




y 




-1 -0.5 





0.5 


1 



. 5 



-0.5 



-1 



Figure 1. Sample eigenvalues of a random matrix 
5^ + VtG w for N = 100 and t = 10" 100 . 

An interesting problem for future research is for a given sequence 
A( N ) f random matrices which converges in ^-moments to x and fixed 
iV to determine the optimal value of for which the measure [Ian) is 

the best approximation of \i x 

Before we present the proofs of these theorems we shall prove the 
following lemma. 

Lemma 8. Let be as in Theorem 0. For every t > and every 
A G C we have that 

lim trln \A[ N) -X\ = In A(x t - A) 
holds almost surely. 

Proof. Let us fix A G C. For any e > we define functions on R + 

ln(r 2 + e) 



Mr) 
g e {r) = In r — 



2 ' 
ln(r 2 + e) 



Each function f e is well defined on [0, oo) and /<= converges to the func- 
tion t i — ► In t pointwise as e tends to 0. Each function g t is increasing 
and g t converges pointwise to as e tends to 0. 
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Figure 2. Sample eigenvalues of a random matrix 
»W + vfa?W for iV = 100 and t = 10" 5 . 
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Figure 3. Sample eigenvalues of a random matrix 
5^ + y/iGW for AT = 100 and t = lO" 2 . 
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Figure 4. Sample eigenvalues of a random matrix 
+ v^G (JV) for JV = 100 and t = 3 • 10 _1 . 

Since the function f e has a polynomial growth at infinity, therefore 
there exist even polynomials S(r) and Q(r) such that S(r) < / e ( r ) — 
Q(r) holds for every r > and furthermore |i?(r) — fe{f)\ < e, |Q(r) — 
/ e (r)| < e hold for every < r < ||xt||. 

We apply Theorem [| for a pair of matrices and AS™ — \ and 
obtain a probability space (O, i3, P) and Gaussian random matrices 
& N \G^ G TWjv^ 00 -^)) such that 

trg £ (\ViGW\)<trg e (\4 N) -\\) 

holds for every a> G n, where as usually A { t N) = + G (JV) . For 
simplicity here and in the following we skip the obvious dependence of 
random variables on uj. 
We have that 

trln \Af > - A| = tr / e (|^ - A|) + tr&fl^ - A|) > 
trSd^ - A|) +tr» e (|Vt = 
tr S(\A[ n) - A|) +trln(|Vt & N) \) - tr f e (\V~t G {N) \) > 
tr S(\Al N } - A|) +trln(|Vt G (iV) |) - trQ(|Vt G w \) =: 1^. 
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Propositions |3| and ^| show that 

Hm XW = <P{S(\x t - A|)) + 0(ln \V~tc\) - 4>{Q{\Vtc\)) > 

0(ln|x t -A|) + 0(^(|Vi c\)) -2e 



holds almost surely. Hence by taking the limit e — > oo we obtain that 
the inequality 

liminftrm|Af° - A| > lnA(x t - A) 

holds almost surely. 
The upper estimate 

lim sup tr In | A[ N ^ — A | < In A (x t — A) 

follows from Lemma and Proposition |3], what finishes the proof. □ 

Proof of Theorem For the proof of the first part of the theorem let 
let K C C be a compact set. In the following y will denote either 
x G A or the matrix A[ N \uo). Let / G C 2 (K) be a smooth enough 
function with a compact support K C C. From the definition of the 
Brown measure we have 

/ /(A) ^(A) = ^(/(A),V 2 lnA( 2/ -A)) = 
Jc 2n 

±- [ \nA(y-X)V 2 f(X) dX 
^ Jc 

Since twice differentiable functions C 2 (K) are dense in the set of all 
continuous functions C(K) therefore the almost certain convergence of 
measures H a (n) ^ in the weak topology to the measure fi x would follow 

if the sequence of functions tr In \A[ N ^ (u) — A| converges to the function 
In A(x — A) in the local C 1 norm almost surely. Therefore it would be 
sufficient to show that for almost every u G Q we have (for simplicity 
here and in the following we skip the obvious dependence of random 
variables on ui) 



lim 

N—KX 



K 



tr In \a\ N ^ - A| — In A(x t — A) 



dX = 0. 



From Lemma [8] and the Fubini theorem follows that for almost every 
wGliwe have 

lim trln \A[ N) - A| = In A(x t - A) 

iV— >oo 
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for almost all A G K. Now it is sufficient to show that 
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(7) 



sup 

K N 



trlii 1 4*° " A l 



dX + 



A 



lnA(art - X) 



dX < oo 



holds almost surely in order to apply the majorized convergence theo- 
rem. 

Note that A i— >• log A(x t — A) is subharmonic ||Bro|| and hence it is a 
local C 1 function; therefore we only need to find estimates for the first 
summand in (0). 

Theorem ^ gives us that for almost every uj G Q 

trm|Af° - A| > trln \VtG^\ 



hence Proposition |9] implies that 



Emin (o, inf trln \Al N) - A|) > Emin (o, inf trln \G {N) \ + — 

is uniformly bounded from below over A G C. From Fubini theorem 
follows that for almost every ueOwe have 



A 



min ( 0, inf tr In \A 



to 



N 



Al )dX > — oo. 



From the simple inequality log r < r 2 which holds for every r > we 
have 



trln \A 



(JV) 



A| < tr \A[ N) - A| 2 < 

W|2 



tr \A 



(JV)|2 



+ A 2 



By Proposition g we have that tr | A\ \ converges almost surely, hence 
the family of functions K 3 X i— > trln — A| is almost surely uni- 
formly bounded from above, what finishes the proof of the first part of 
the theorem. 

From the first part of theorem follows that there exists a decreasing 
sequence (%) of positive numbers which converges to and such that 
for any compact K G C 



lim 



K 



trln | A™ 



X\-\nA(x tN - X) 



dX = 



holds almost surely. 

Proposition ^ implies that the majorized convergence theorem can be 
applied (we recall that A i— > log A(y — A) is always a local C 1 function) 
hence 



lim 

N—KX 



K 



In A(x tN - A) - In A(x - A) 



dX = 0. 



The above two equations combine to give 



lim 

TV^oo 



A' 



trln|Ag° - A| - In A(x - A) 



dX = 
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almost surely. The convergence of empirical distributions of eigenvalues 
follows now exactly as in the proof of the first part. □ 

Proof of Theorem Let (tjv) be a sequence given by Theorem |. Since 
limsupjv^oo HG^H < oo almost surely ||Gem|| , hence for almost every 

is the wanted sequence. □ 

5. Technical results 
5.1. Derivation of the stochastic differential equation for A«. 

5.1.1. Singular values as functions on M. N . In this subsection we are 
going to evaluate the first and the second derivative of the map 

s : Mn ^ (si(m), . . . , SAr(m)) , 

where Si(m), . . . , ^(m) denote the singular values of a matrix m. 
The perturbation theory shows (cf chapter II. 2 of |Kat| ) that if D is 



a diagonal matrix with eigenvalues v\ , . . . , such that ^ Uj for all 
i 7^ j and AD is any matrix then the eigenvalues u[, . . . , v' N of a matrix 
D + AD are given by 

v\ = u t + AD U + AAj _ AZ ^ + 0{\\ADf) 

... i^i Vj 

3 A 

for small enough || AF|| and that the map AD \— > . . . , v' N ) is C 2 in 
some neighbourhood of 0. 

It follows that if F is a diagonal matrix with positive eigenvalues 
Si,...,Sjv, and AF G M.n is any matrix then the singular values 
s[, . . . , s' N of F + AF are given by 

2 



(8) (s' l ) 2 = s* + 2sMAFu + J2\ AF ; 



3 

3+i ~ ^ 

and that the map AF i— > (s[, . . . , s' N ) is C 2 on some neighbourhood of 
0. 

In the general case every matrix X can be written as X = UFV, 
where F is a positive diagonal matrix and U, V are unitaries. If the 
singular values of X are s\, . . . ,sn then (||) gives us singular values 
s[, . . . ,s' N of the matrix X + AX, where AF is defined by AF = 
V*AXU*. 
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5.1.2. Trajectories of the Brownian motion avoid sigularities of s. The 
set of singularities of the map s, namely 

{m G M. jv : det m = or Sj(m) = Sj{m) for some i ^ j}, 

is a manifold of codimension 2 and hence almost every trajectory of a 
matrix Brownian motion A t will avoid this set. In this subsection we 
will present a rigorous proof of this statement. 
For every e > we define a set 

K e = < m G M.n '■ ^^lnsi(m) > e, ^^hi \si(m) 2 — Sj(m) 2 \ >e>. 

^ i i<j * 

First of all, for any fixed e > we define a stopping time 
T(u) = min{t > : In | det A t (u) \ < e} 
and a stopped Brownian motion 

A t {u) = Anin[t,T(a;)](^)- 

In the proof of Proposition [| we showed the function m t— > In | detm| 
is subharmonic and hence 1 1— > In | det A t \ is a submartingale. It follows 
that 

(9) In | det A\ < Eln|deti T | < Emax (0, In | det A T \) + 

lne P(uj G n : In | det A t (oo)\ < e for some < t < T). 

Since the first summand on the right-hand side of the above inequality 
is clearly finite, it follows that 

lim P(u GO: I det A t (u)\ > e for every < t < T) = 1. 
Secondly, we consider a function on Mtv(C) given by 

(10) m i — > In | (gj(m)) 2 — (sj(m))' 

Formula (§) gives us first and second derivatives of the map m i— > 
(si(m), . . . , sjv(w)) and allows us to find the Laplacian of the each 
summand in (|T0|): 



V 2 m|^-s?| sj + s) sJ + 4 s) + s\ 

- s 2 ) 2+ ^ {si - s 2 )(s 2 - si) (si - s)){s* - si) 
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It is not difficult to see that for every i,j, k all different we have 



4 + 4 



4)(4 



4M 




4 + 4 



4 + 4 



s .2 



4 + 4 



(4-4)(4-4)J \(4-4)(4-4) (4-4)(4~4). 

and due to these cancellations 



v 2 v>| 3 ?- s ;i=4y; 



i<j 



i<j 



4 + 4 
4 - *?) a 



> o 



holds. It follows that 

t^^ln|( Sl (A)) 2 -(^(A)) : 



i<j 



is a submartingale and by similar arguments as in (|9|) we see that 
lim P(lu £ ^ : y^ln |g 



iym) 2 — Sj(m) 2 \ > e for every < t < T) — 1. 



5.1.3. Stochastic differential equation for Aj. We recall that every ma- 
trix m can be written as m = U(m)F(m)V(m), where F(m) is a 
positive diagonal matrix and U(m), V(m) are unitaries. Let us define 
now a new matrix-valued stochastic process B given by a stochastic 
differential equation dB = (V (A t ))* {dM)U (A t )* . It is easy to see that 
B is again a standard matrix Brownian motion. 

For any fixed e > let us consider any C 2 function s : A4n — » 
R such that s(m) = s(m) for every matrix m £ K e and such that 
|| 5(777.) || < C||m|| for some universal constant C and all m G M.n, 
where || • || denotes any norm on A4 n or M. N respectively. 

Function s = (s 1} . . . , Sjv) fulfills assumptions of Ito theorem, hence 
for every T > we can write the Ito formula 



Si (At) 



E 



dsi{m) 



dm 



ki 



dM kl + 



E 



m=A t 

1 / d 



'2 7. 



Si(m 



AN \(d^m k iy 



+ 



d 2 sAm) 



m=A t 



m=A t 



dt. 



For every u G Q such that A t (u) G K e for all < t < T the left-hand 
side of this equation is equal to Aj(T) and the right-hand side can be 
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computed from 



2A~ 



1 - 



2N + ^ 



Since almost every oj has the property that for some e > we have 
At G K e for all < t < T hence the above equation holds without any 
restrictions for oj. Equivalently, 



d\{t) = M(dB ti ) + 



dt 
2A~ 



1 



2N + ^ 

3& 



A^ + A| 



5.2. Determinant of a standard Gaussian random matrix. 

Proposition 9. Let (G^) be a sequence of independent standard 
Gaussian random matrices and let c be a circular element. Then 



lim trln|C7 (Af) | 

N->oo 



lnA(c) = — 
v ' 2 



holds almost surely. 

Furthermore for any s G 



we have 
Emin (s,inf trln |C7 (A °| 



A' 



> 



-oo. 



Proof. The square of a circular element is a free Poisson element with 
parameter 1. The probability density of this element can be explicitly 
calculated ||VDN|| and the integral lnA(c) = J °°lnr dfi^—^ can be 
computed directly. 

Let us fix iV G N. Let v\, . . . , Vn be random vectors in which are 
defined to be columns of the matrix G^ N '. We define 



Vi 



\ 



det 



(vi,vx) 



(Vi,Vi) 



(Ui.Ul) ••• (Vi,Vi) 

where (■,•) is the standard hermitian form on C N . The above matrix 
[(vk, Vi))i<k,i<i is the complex analogue of the Gram matrix; therefore— 
informally speaking — we can regard Vi to be the "complex volume" of 
the "complex parallelepiped" defined by vectors v\, . . . , 

Of course V i+ i is equal to the product of Vi and l i+ i, where l i+ i 
is the length of the projection of the vector onto the orthogonal 
complement of the vectors v±, . . . , v j. Since 

\Vl,Vx) ■■■ (vi,v N ) 



. (vn,vx) 



(vn,v n ) 



{G (N)y G (N) 
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it follows that 

| det G {N) \ = V N = hh---l N . 

It is easy to see that the distribution of /j coincides with the distri- 
bution of the length of a random Gaussian vector with an appropriate 
covariance in the complex (N — i + l)-dimensional space and therefore 

-h Jo°° r ~ V 2 ( JV - i+1 )- 1 e- JVr2 dr N^T(N - i + 1 - §) 
EL — 



J °° r 2(7V- l+ l)-l e -iVr2 dr ' r (AT _ i + 1) 

and hence Markov inequality gives us 

P [(h ■ ■ ■ In-i)- 2 > e (1+2£)7V ] < e (-i-*)"tf"-i * < e - 2 ^. 

Above we have used that random variables ij are independent and 
simple inequality r(JV) > (^P 1 )^ 1 for iVeN. 
Since 

p( 1 ° g/l + - A ; + lQg ^<-l-2e)< 

/logZi + --- + logZ JV _i 1 \ /log /at ^ 

Borel-Cantelli lemma implies 

liminftrln|G (7V) | > lnA(c) 

almost surely. This together with Lemma ^ gives us the first part of 
the proposition. 

It is possible to find a constant C such that for every e > ■? we have 
P ^inftrlnlG^I < ~ - 2eJ < ^ F ^tr In | G (iV) | < -| - 2eJ < Ce" 

If z/ is the distribution of the random variable inf jv tr In | | then 
integration by parts gives 

/ (t + l)dv{t) = - j p(-oo,t) dt > -oo 

J — oo J~ oo 

and the second part follows. □ 
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